Controlled manipulation of light by cooperative response of atoms in an optical lattice 
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We show that a cooperative atom response in an optical lattice to resonant incident light can be 
employed for precise control and manipulation of light on a subwavelength scale. Specific collective 
excitation modes of the system that result from strong light-mediated dipole-dipole interactions can 
be addressed by tailoring the spatial phase-profile of the incident light. We demonstrate how the 
collective response can be used to produce optical excitations at well-isolated sites on the lattice. 
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Accurate control of ultracold atomic gases in periodic 
optical lattices, in which interactions are well understood, 
opens the door to unique and intriguing opportunities to 
study many-particle phenomena and their applications. 
Experimental progress has lead to observations of novel 
strongly-interacting states, e.g., in quantum phase tran- 
sitions [H-Q a-nd fermionic pair condensation Q. Many- 
body quantum entanglement has been generated via con- 
trolled atom collisions , lattice systems have been used 
for preparation of spin-squeezed states for sub-shot-noise 
interferometry Q, and the atoms can now even be ma- 
nipulated in a single-spin level at a specific lattice site 
[9|. On the other hand, recent developments in nanofab- 
rication of arrays of circuit elements acting as plasmonic 
resonators has stimulated interest in photonic metama- 
terials. A metamaterial is an artificially tailored crystal 
consisting of subwavelength-scale structures that can ma- 
nipulate light on a nanoscale. Here we show that a basic 
Mott-insulator state of a neutral gas of ultracold atoms 
confined in an optical lattice, or artificial light crystal, 
exhibits strongly interacting electric dipole transitions 
leading to a cooperative response. Such collective be- 
havior can influence resonant imaging and may also be 
employed to form a metamaterial for precise control and 
manipulation of optical fields on a subwavelength scale, 
providing an interesting nanophotonic tool. 

In this letter, we consider an ultracold gas of atoms 
confined in a two-dimensional (2D) optical lattice with 
precisely one atom per site. Such a system can be pre- 
pared, e.g., in a weak harmonic trap or by engineer- 
ing a Mott-insulator state of atoms by single-site ad- 
dressing of atomic spins Q. Resonant, coherent light 
whose spatial phase-profile is adjusted, e.g., by a holo- 
gram or spacial light modulator, illuminates the lattice. 
The scattered light mediates strong many-particle dipole- 
dipole (DD) interactions between atoms, leading to a 
cooperative atom response. The optical excitations of 
the atoms exhibit collective modes with resonance fre- 
quencies and radiative linewidths that dramatically dif- 
fer from those of an isolated atom. We demonstrate the 
idea of subwavelength-scale light manipulation by engi- 
neering the spatial phase-profile of an incident monochro- 
matic plane-wave. The tailored incident field produces 
localized dipolar subwavelength-scale excitations of the 
atoms in desired locations in the lattice. By dynamically 



adjusting the phase pattern of the field, the excitations 
can be controlled and moved around in the lattice. 

The particular example of subwavelength-scale local- 
ization of optical excitations we study here has attracted 
considerable interest in nanophotonics with possibili- 
ties for microscopy and data storage applications. Spa- 
tial and temporal modulation of ultrashort laser pulses 
leads to excitation of energy hot-spots in nanostructures 
[lOl [Tlj . It has also been proposed that interactions be- 
tween induced currents and plasmonic waves on nanos- 
tructures permits the excitation of subwavelength hot- 
spots by amplitude or phase modulated monochromatic 
fields [lllil. 

Non-trivial collective optical properties result from a 
cooperative response of the strongly interacting, closely- 
spaced atoms: recurrent scattering events, in which a 
photon is repeatedly scattered by the same atom, lead 
to collective modes with strongly modified spatial con- 
figurations and radiation rates p4l - [l7j . Such scattering 
processes can result in light localization that is analogous 
to Anderson localization of electrons [l^. The resonant 
response is very different from the studies of off-resonant 
optical diag nostics of atomic correlations in optical lat- 
tices |19l - [25| . Photonic band gaps for atomic lattices have 
previously been calculated in Ref. [1^ . 

For our lattice system we numerically calculate the op- 
tical response by stochastically sampling the atomic po- 
sitions according to their spatial distributions then solv- 
ing the recurrent scattering events to all orders for each 
stochastic realization. We find a strong resonant response 
in the case of closely-spaced atoms with the near-field 
emission pattern from the atoms forming sharp, nar- 
row spatially localized amplitude peaks. The response 
is sensitive to detuning of the incident light from the 
atomic resonance and to the spatial separation between 
the atoms. Tuning light off-resonant or increasing the 
lattice spacing rapidly leads to suppressed interactions. 
If the atoms are not confined strongly enough to the in- 
dividual lattice sites, the resulting increased disorder in 
the atomic positions due to quantum fiuctuations also 
suppresses the strong collective effects in the ensemble- 
averaged response. 

We take the atoms to occupy the lowest energy band in 
a 2D square optical lattice of periodicity a in the xy plane 
as described in Appendix \^ We assume that the atoms 
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are tightly confined to the lowest vibrational state in the 
z direction of an oblate external potential and that they 
reside in the Mott-insulator state with precisely one atom 
per site. In a combined harmonic trap and the lattice the 
single-occupancy state of bosonic atoms can exist in a 
weak harmonic trap or can be engineered from the typical 
'wedding-cake' Mott-insulator ground state by manipu- 
lating the multi-occupancy states, e.g., by single-site ad- 
dressing or by inducing atom parity-dependent losses . 
In a deep lattice, the vibrational ground-state wavefunc- 
tion (the Wannier function) is approximately that of a 
harmonic oscillator with frequency w = where 
s denotes the lattice depth in the units of the lattice- 
photon recoil energy = ■k'^ / {2ma?) [l^]. The con- 
finement of the atoms along the lattice in the xy plane 
in each site j thus has a Gaussian density profile, Pj{r), 
with the 1/e width (. = as^^l'^ jix controlled by the lattice 
depth s and the lattice spacing a. 

We illuminate the lattice with a monochromatic in- 
cident field Ein(r, whose frequency 57 is nearly res- 
onant on an electric dipole transition. This impinging 
field excites the dipole transition of the atoms, produc- 
ing scattered light that, in turn, impact the driving of 
neighboring atoms and alter their scattered light. The 
scattered photons can mediate strong interactions be- 
tween closely-spaced atoms, so that the atomic system 
responds to light cooperatively, exhibiting collective ex- 
citation eigenmodes. Here, we show how to exploit these 
interactions for controlling and manipulating light on a 
subwavelength scale. As a specific example, we prepare 
subwavelength-scale spatially localized collective excita- 
tions of the atoms in isolated regions of the lattice by con- 
sidering an incident plane- wave illumination of the atoms 
with an approximately sinusoidal phase-profile. Such a 
response is distinct from that which would be seen if the 
atoms did not interact. 

In order to model the cooperative atom response to 
light, we assume the incident field is sufficiently weak 
that saturation of the excited state can be neglected. For 
simplicity, we consider the atomic internal states as an 
effective two-level system consisting of a single electronic 
ground and excited state. On impact, light drives atomic 
transitions inducing a polarization density P"'"(r, — 
Pj*" (r, i) , where the polarization within each site j, 
P^(r,t) = dj(5(r — rj) and dj is the electric dipole mo- 
ment of an atom at site j. To facilitate numerical evalu- 
ation of the lattice response, we express the polarization 
in terms of stochastic amplitudes Cj , representing the co- 
herence of atoms 2 realized for a stochastic sampling of 
atomic positions from the atomic density distributions 
Pj{rj), such that P^(r, = e~'^*pdJ(r—rj)ej(t), where 
p is the atomic dipole matrix clement. When the atomic 
dynamics evolve on timcscalcs much longer than the light 
propagation time across the optical lattice [l6j . the in- 
duced polarization produces the scattered electric dis- 
placement field D+^. (r) = fcV(47r) / d3r'G(r-r')-P+(r'), 
where G(r — r') is the monochromatic dipole radiation 
kernel representing the radiated field at r from a dipole 



residing at r' [28[. Thus, the atom at site j experiences 
driving by the sum of fields Dg j' scattered from all other 
atoms in the lattice and the incident field Din; these 
scattered fields are proportional to the amplitudes eji of 
their atoms of origin. These multiple scattering processes 
therefore produce collective dynamics described by 

e, = {i6 - r/2)e, + ^ C,,,,e,, + F„ (1) 

C„v ^ — d . G(r, - r,0 • d, (2) 

where Fj = e^^*' pd-'D^^{rj) / {iheo) is the direct driving of 
the atom in site j by the incident field, S = fl — uje.g is the 
detuning of the field from resonance, and T is the spon- 
taneous emission rate. For each stochastic realization of 
atomic positions, interactions between N atoms in an op- 
tical lattice lead to the formation of N collective atomic 
excitations, each with its own resonance frequency and 
spontaneous emission rate, which could have either su- 
perradiant or subradiant characteristics. 

By exploiting the strong DD interactions, one can tai- 
lor the incident field so that it excites specific linear 
combinations of collective modes, providing a desired 
response. Here, for example, we consider a linearly- 
polarized {ey chosen to be parallel to the lattice plane), 
phase-modulated field whose positive frequency compo- 
nent reads 

E+(r,i) ~ eyi;oe"^(^''')e*('^'^-^"\ (3a) 

TT 

(p{x, y) = — sin(Ka;) sin(K?/) . (3b) 

A field profile of this form, however, contains evanescent 
plane- wave components whose transverse wavevectors ex- 
ceed the carrier wave number k = 27r/A = fi/c. We 
therefore approximate it through the truncated Fourier 
expansion as described in Appendix [Cl One can produce 
such a field, e.g., by a hologram or a spatial light modu- 
lator. Figure [T] illustrates how the field Eq. ([3]) can excite 
a checkerboard pattern of localized excitations of atoms 
in a lattice. In these calculations, we neglect the width of 
the atomic wavefunctions in the z direction, as a width 
on the order of I in the z direction has a negligible ef- 
fect on the calculated response. We fully incorporate, 
however, quantum fiuctuations of the atomic positions 
on the lattice plane in the vibrational ground state of 
each site. The many-atom correlations of the specific 
one-atom Mott-state are also included. As we discuss 
below, interactions between the atoms are vital to the 
realization of this pattern. The atoms are arranged in an 
18 X 18 square lattice with the spacing a = 0.55A, and a 
site residing at (a;o,yo) = (a/2, a/2). The incident field 
[Eq. Q] has a modulation period 2t: / k = 6a, indicating 
the periodicity of Ein(r,i) of six sites. We choose the 
dipole orientation A k, + O.le^ to be slightly rotated 
from the normal to the lattice plane so that the atoms 
scatter fields largely within the plane while allowing them 
to be driven by the incident field. The collective mode 
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FIG. 1: The atomic excitation intensities \e.j\^ I maxj/ |ej/P) resulting from the cooperative response of an ensemble of two-level 
atoms in an optical lattice to a phase modulated incident field [Eq. ((3)l] with modulation period of 6 lattice sites. The black dots 
indicate the positions of the populated lattice sites. In panel (a), the atoms are perfectly confined at the center of each site and 
their wavefunctions have width £ = 0. Panel (b) represents the response of a single stochastic realization of atomic positions 
sampled from a 2D Gaussian variables of width 0.12 a centered on the lattice sites. Panels (c), (e-f) the show responses of 
atoms whose wavefunction widths £ = 0.12 a calculated from an ensemble average over several thousand realizations of atomic 
positions, while in panel (d) £ = 0.21a. In panel (e), the incident field is detuned from the resonance of an isolated atom by 
lOr, while the detuning is zero for all other panels. In panel (f), atoms were removed from 8/9 of the sites, resulting in an 
effective spacing of a' = 3o, thus removing the role of the other sites from the cooperative response and destroying the pattern 
observed in panels (a) and (c). 



spontaneous emission rates range from a very subradiant 
3 X lO^'^r to the superradiant 5T, while their frequency 
shifts from the single atom resonance range from — 2r to 

o.9r. 

We first consider an infinitely deep lattice in which the 
atoms are perfectly confined at the center in their re- 
spective sites, i.e., with i = 0. In this case, the atomic 
positions are deterministic and Eq. ^ reduces to a cou- 
pled set of linear equations whose steady-state solutions 
for |ej p are shown in Fig.[T^. A subset of atoms residing 
at the local minima of ip(x,y) are more strongly excited 
than those in the surrounding lattice sites, while those at 
the local maxima of (p{x, y) are roughly as weakly excited 
as their surroundings. The peaks sit on a background 
with excitations roughly 0.2 times those of the most ex- 
cited atoms. We find a subwavelength excitation FWHM 
width of the peak to be less than 0.9A. This results in a 
checkerboard pattern with the localized excitations sepa- 
rated by 3\/2 sites (3 sites in the both directions) sitting 
on a background of weakly excited atoms. The regular- 
ity of the response can be broken by slightly altering the 
period or the orientation of the phase modulation. Note 
that the periodicity of the incident field is significantly 
larger than the width of the localized excitations. 



In a more realistic scenario, the width of the lattice site 
atomic wavefunction due to vacuum fiuctuations cannot 
be neglected and the scattered light sources are essen- 
tially distributed over the atomic densities. To obtain 
the average atomic excitations |e|| that dominate the 
near-field emission, we solve Eq. ([1]) through Monte-Carlo 
integration [2^. We obtain a large number of realiza- 
tions of atomic positions Vj in each site j sampled from 
a probability distribution matching the single-atom den- 
sity function. Then for each realization, we solve Eq. ([1]) 
as if the atoms were localized at the sample points. We 
then compute |ejp for each sample, and perform an en- 
semble average over all realizations of atomic positions. 
Each stochastic realization represents a possible outcome 
of a single experimental run in which case atomic posi- 
tions are localized due to the measurement of scattered 
photons. 

Quantum fiuctuations of atomic positions in individual 
sites can dramatically affect the response. We demon- 
strate this showing in Fig. [T]3 the excitation intensities 
for a single stochastic realization of atomic positions in 
which case a single atom in each site is independently 
sampled from the Gaussian density-distribution of width 
i = 0.12 a (corresponding to the lattice height s ~ 50). 
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The non-regularity of the lattice alters the collective 
modes for the sample realization. Although the atoms 
are driven by the same incident field that gives rise to 
the pattern in Fig.[T^, the displacement drastically alters 
the collective interaction, and yields an optical response 
with a significant stochastic noise and a less regular array 
of localized peaks. Such effects can be washed out when 
one calculates the ensemble average of the response that 
corresponds to expectation values obtained over many ex- 
perimental realizations. The excitation intensity je^P av- 
eraged over 6400 position realizations for Wannier func- 
tions of width I = 0.12 a is shown in Fig. [TJ: for the 
same parameters as those in Fig. [1^. With atomic wave- 
function of this width, the collective interactions produc- 
ing the pattern of Fig. [1^ survive the averaging process, 
providing an excitation with subwavclcngth FWHM and 
a background excitation comparable to that calculated 
for perfectly localized atoms in an infinitely deep lattice. 
Weaker confinement only moderately diminishes the co- 
operative interactions. For i = 0.21a, corresponding to 
s ~ 5, Fig. [TJl shows a weakening contrast of the pattern 
with a background excitation approximately 0.3 times 
that of the peaks, which themselves have slightly broader 
FWHM of 1.2A. 

We can illustrate the essential nature of cooperative in- 
teractions in the formation of this excitation pattern by 
suppressing the DD interactions and observing the result. 
Fig. [TJ) shows the response of the atoms in the optical 
lattice of Fig. [TJ; to an incident field detuned from the 
single atom resonance by lOF. Since, with this detuning, 
each photon interacts more weakly with an individual 
atom, the probability of a multiple scattering process for 
a two-level atom reduces geometrically with the number 
of single scattering events involved in that process p7| . 
The interactions resulting from these multiple scattering 
events are therefore suppressed. While the incident field 
has the same phase modulation as in Figs. [TJi-d, the ex- 
citation pattern is not preserved. Atoms at sites where 
the phase modulation ip is minimized are not apprecia- 
bly more excited than the surrounding background as in 
Figs. [1^ and c. 

The pattern formation in Figs. [TJi,c truly represents 
a cooperative response where the interactions between 
all the atoms in the lattice, including also the weakly 
excited ones, are essential. We illustrate this in Fig. [If, 
where all of the atoms not residing at the maxima or min- 
ima of the incident field phase modulation ip have been 
removed from the lattice. In effect, the lattice spacing 
was tripled to a' = 3a = 1.65 A, with an atom resid- 
ing at (xo,yo) = (a'/^iO'/S). If the removed atoms had 
not played an essential role, the response of the lattice 
would show a checkerboard pattern of strongly excited 
and weakly excited atoms. However, the system response 
displayed in Fig. [If shows that the atoms in the interior 
of the lattice are excited roughly evenly even though each 
atom is driven with an opposite phase to that of its near- 
est neighbor. Atoms at the edge of the sample are more 
strongly excited due to finite size effects. 



Sharp localized excitations may be broadened by heat- 
ing and losses that can inhibit the cooperative atom re- 
sponse. Raman transitions to other vibrational center- 
of-mass states heat up the atoms, broadening the atomic 
density distributions in individual sites and increasing 
the hopping amplitude of the atoms between the adja- 
cent sites. Such processes could be reduced, e.g., due to 
the orthogonality of the cigenfunctions in each site, if the 
electronic ground and excited state atoms approximately 
experience the same lattice potential even if the system 
is not in the Lamb-Dicke regime. Alternatively, if the 
atomic linewidth is much larger than the trapping fre- 
quency, the collective response may reach a steady state 
before the heating becomes deleterious. Additional atom 
losses due to Raman transitions to other internal states 
could be avoided if each atom forms an effective two-level 
system of a single electronic ground and excited state. A 
desired two-level configuration could be realized with a 
cycling transition by shifting all other transitions out of 
resonance or with aJ = 0— ^J^l transition with all 
except one excited state shifted out of resonance. 

In conclusion, we have shown that resonant DD in- 
teractions between atoms in an optical lattice lead to a 
collective response that can be exploited in manipulation 
of light on a subwavelength scale. To illustrate this, we 
studied an example of engineering a checkerboard pat- 
tern of isolated atomic excitations. Unlike in nanofab- 
ricated metamaterial samples (lol - [l3| . here the effect is 
not based on interactions between plasmonic and current 
excitations but purely electric DD interactions between 
neutral atoms without magnetic contributions. More- 
over, the positions of the excitations can be dynamically 
altered simply by translating the phase-modulation pat- 
tern, so that the collective excitation pattern adiabati- 
cally follows the change in the phase pattern. By un- 
derstanding these interactions, the characteristics of an 
incident field could be engineered to produce more com- 
plex excitations. Our example also demonstrates how a 
cooperative response can have implications on the reso- 
nant absorption imaging of 2D atomic samples in which 
case deviations from the column density results have been 
experimentally observed [3Ci |. Moreover, the non-trivial 
relationship between incident field modes and the col- 
lective excitations to which they couple could be of im- 
portance, e.g., to imaging and to the implementation of 
quantum memories in optical lattice systems in which the 
states of a light field are mapped onto collective hyperfine 
excitations in the lattice. 
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Appendix 

In these appendices, we provide some of the technical 
details of the optical lattice system and the incident field 
used to excite it. In Appendix [A} we describe the optical 
lattice potential. In Appendix[Bl we describe the dynam- 
ics of the atomic dipoles interacting with both the inci- 
dent and scattered electric fields. Finally in Appendix [C| 
we describe how one can construct the approximate phase 
modulated driving field used to illuminate the optical lat- 
tice system and excite collective atomic excitations that 
were discussed in the text. 



Appendix B: Dynamics of atomic excitations in the 
lattice 

In this section, we elaborate on the dynamics of the 
atomic electric dipoles interacting with light. We show 
that the evolution of the system can be described as in 
Eq. (2) of the main text. We then discuss the Monte- 
Carlo m ethod used to compute the average excitation 
energies |ejP displayed in Fig. 1 of the text. 

1. Basic model 



Appendix A: The Optical Lattice Trapping Potential 

We consider a two-dimensional square optical lattice 
of periodicity a, centered on the point {xo,yo,0). Four 
intersecting beams produce the optical confining poten- 
tial 



V = sEr 



2 1 X- Xq 
sm TT 



2 1 y-yo 

sm TT 



We consider an ensemble N two-level atoms placed in 
harmonic potentials centered at the lattice sites Rj (j = 
1, . . . ,N). Each atom has an internal ground state \g) 
and an excited state |e) which differ in energy by huJe,g- 
A coherent, monochromatic field Dj^(r, t) with frequency 
VL impinges on the lattice. The field incident on the atoms 
in the lattice drives the \g) o |e) transition, inducing a 
polarization density P+(r,t) = X]j where the 

(Al) polarization due to an atom in site j, P^^ = dj(5(r — r^), 
dj is the dipole moment of an atom in site j, and r. 



where Er = tt'^U^ / (2TOa^) is the lattice recoil energy pTf . 
and the dimensionless confinement strength s controls 
the width of the vibrational ground-state wavefunction 
in each site (Wannier functions) that results from single- 
particle quantum fluctuations. An additional potential 
tightly confines the atoms in the 2 = plane. The lat- 
tice resides in a Mott-insulator state with precisely one 
atom per lattice site. Each lattice site, labeled by in- 
dex j , has a potential minimum located at position Rj, 



and the Wannier function 4>j{Y) 



(r - Rj). When 



the confinement is sufficiently tight, the potential V in 
the neighborhood of each lattice site is roughly harmonic 
with 



(A2) 



where 



2^Eii/h, uJz is the trapping frequency 



along the z direction, and Ar^ = • (r — Rj) is the 
displacement of the //th component of r from the lattice 
site Rj . An atom in each site occupies the ground state 
of the harmonic oscillator potential 



(r) 



1 



(^3£4^2)l/4 



exp 



y 



2P 



2P 



(A3) 



where the width of the of the wave function is £ = 
as~^/'^ /tt, and its thickness = \/ h/{mujz). The atomic 
density, Pj{v) = |0j(r)p, at site j thus has a Gaussian 
profile with a 1/e radius i in the xy plane. This width 
is directly proportional to the lattice spacing and nar- 
rows with increased trapping strength s. It is possible 
to adjust the lattice spacing by using accordion lattices 
jsif . Moreover, the relationship between the lattice spac- 
ing and the Wannier wavefunction confinement may be 
controlled by magnetic field dressing [32j . 



is its position coordinate. Scattering of light from the 
atoms produces correlations between the positions of the 
atoms and the scattered field profile. In this way, contin- 
uous measurement of the scattered light yields an effec- 
tive measurement of the atomic positions. This process 
simultaneously projects the position coordinates of an 
atom in each site j onto a specific positon through 
the projection operator Q{{vj}) = (S'jLi acting 
on the initial vibrational state of the system, where {tj } 
denotes the set {ri, . . . , tat} of realized position coordi- 
nates for atoms in lattice sites {1, . . . , N}. We therefore 
characterize the evolution of the polarization density as- 
sociated with site j in terms of the slowly varying coher- 
ence conditioned on the realized atomic position Vj , and 
the positions in all other sites Yji for j' ^ j. We define 
the amplitudes ej{rj,t]{vji^j}) = e'^^*((Tj){ri,r2,...,riv} to 
be the slowly varying coherence between the ground and 
excited states of an atom at site j conditioned on the 
observation of atoms at positions {ri, r2, . . . , r^} within 
their respective sites, where Uj is the coherence operator 
for an atom in site j. For a single realization of position 
coordinates of the atoms, the polarization in site j de- 
pends on the position of all the atoms via the amplitudes 
ef 

P,(r,t;{r,v}) = e-'"*dp5(r-r,)e,(r„ t; {r,v,}), (Bl) 

where p is the dipole matrix element. In addition to 
the incident electric field, an atom at site j experiences 
driving by the fields Dg scattered from all other sites 
f 7^ j in the lattice. The total displacement field arriving 
at this site is thus 



D 



GXt ,j 



-D+(r,t) + ^D+^.,(r,t). 



(B2) 
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When the atomic dynamics evolve on timescales much 
longer than the time it takes for light to propagate across 
the optical lattice, the scattered fields can be expressed 
in the monochromatic limit as [l7| 



Dj,,(r) = |^ I dVG(r-r')-P+(r'), 



(B3) 



where k = U,/c, and G is the radiation kernel with tensor 
components 



1 



(r) = p {d,d, - ^M.^V^) 



(B4) 



The monochromatic dipolc radiation kernel, representing 
the radiated field at r from a dipole residing at r', can 
be expressed as [11] 



G(r) 



2 /rr 1 

3 ^ 3 
Stt 



3 - 3ikr - k'^r'^ 



^ikr 



-5{kv). 
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Stochastic simulations 



In the simulations we solve the cooperative optical re- 
sponse of the atomic sample to the incident field using a 
Monte-Carlo approach in which the position coordinates 
of the atoms are sampled according to their position cor- 
relation functions, and the optical response is calculated 
for each stochastic realization [2^. We assume that the 
atoms form a Mott-insulator state of precisely one atom 
per lattice site. The position coordinates of atoms within 
each site j are therefore independent stochastic vari- 
ables, sampled from Gaussian distributions Pj(r). The 
atoms reside at the vibrational ground states of the lat- 
tice sites and the Gaussian distribution pj (r) results from 
the single-particle quantum fluctuations, determined by 
the Wannier wavefunction density in the particular site j . 
The atomic position in the site j is centered at Rj . The 
width of the Wannier wavefunction is determined by the 
lattice confinement. Since each stochastic realization of 
position coordinates for the atoms in the N lattice sites 
{ri,r2, . . . ,rj^} can be interpreted as an outcome of a 
continuous measurement process of scattered light that 
localizes the atoms, each stochastic Monte-Carlo trajec- 
tory also represents a possible outcome of a single ex- 
perimental run. We evaluate the ensemble averages of 
the atomic optical excitations by computing the excita- 
tions for each of many stochastic realizations of position 
coordinates of the atoms and then calculating their en- 
semble average. The numerically calculated ensemble av- 
erage then corresponds to the experimentally measured 
ensemble average over many experimental runs and pro- 
vides quantum mechanical expectation values of the ob- 
servables. 

In order to calculate the optical response for each 
stochastic realization of position coordinates for the 



atoms in the N lattice sites, {ri, r2, . . . , r^v}, we assume 
here that the incident light is sufficiently weak that we 
can neglect the excited state saturation of the atoms. In 
this limit, the coherence amplitude for an atom in site j 
obeys the equations of motion 



-e, = {i6 - T/2)e, 



s'"*-£-d.D+,^^.(r„t), (B6) 



where S = fl 



is the detuning of the incident light 



from the resonance of the atomic transition and F is the 
atomic spontaneous emission rate. For each stochastic 
realization of atomic positions, the light impinging on an 
atom in a particular site consists of the incident field and 
the scattered light from all the other iV — 1 sites. Because 
light scattered from an atom is directly proportional to 
its coherence amplitude, we obtain 

= (z5-r/2)e, + ^C„.e,.-fe^^^*-£-d.D+(r„t), 



(B7) 

where Cjj' is the coupling matrix element representing 
the effect of dipolc radiation from the atom in site j' on 
the coherence amplitude of the atom in site j ; 



■ G(r, 

2i ^ ^ 



r,0 . d. 



(B8) 



For each stochastic realization of position coordinates of 
the atoms, Eq. (jB7p represents a collective response of the 
atomic sample to the incident light. Repeated exchanges 
of a photon between the same atoms lead to coopera- 
tive response of the atoms. The system exhibits collec- 
tive eigenmodes with each eigenmode associated with a 
specific resonance frequency and a radiative lincwidth. 
In the text, we considered an 18 x 18 array with lat- 
tice spacing 0.55A (A — 27rc/r2), and dipole orientations 
d « -I- O.ley. For this configuration, when the atoms 
are placed precisely at the centers of the lattice sites 
Hj, the collective mode spontaneous emission rates range 
from a very subradiant 3 x 10~'^r to the superradiant 
sr. while their collective frequency shifts from the single 
atom resonance range from — 2r to 0.9F. 

In each stochastic realization of position coordinates 
of atoms within their respective lattice sites, we solve 
the steady-state solution of the excitation amplitudes of 
the atom in each lattice site Cj by setting dcj/dt = in 
Eq. (|B7p . The results for a single stochastic realization 



of \ej\ are displayed in the main text. We then compute 



the ensemble average of the atomic excitation 



by 



averaging over a sufficiently large number of stochastic 
realizations. That is, we take Af stochastic realizations of 
the position coordinate of an atom within each site j. For 
the i/th realization (v = 1, . . . ,Af) we have sampled the 
position coordinates of an atom in site j, r^"''. For these 
atomic positions, we solve the steady state of Eq. (jB7[) 
to obtain the atomic amplitudes e^'''' and calculate the 
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ensemble average 

v—l 



3. Many-particle correlations 



In the resulting joint probability distribution the posi- 
tion coordinates of atoms within each site j are indepen- 
dent stochastic variables that are sampled from the Gaus- 
sian distributions Pj(r) = |0j(r)p, corresponding to the 
stochastic Montc-Carlo sampling procedure implemented 
in the previous section. 

Appendix C: The incident light 



In this paper we consider a specific many-particle atom 
state: the bosonic Mott-insulator state with precisely one 
atom per site. The atoms reside in the lowest vibrational 
state of each lattice site. Vacuum fluctuations for the 
position of the single-particle state of the atoms are in- 
corporated in the Monte-Carlo sampling of the atomic 
positions from the Gaussian density distribution, deter- 
mined by the Wannier wavcfunctions of each lattice site. 
As we will show below, for the specific Mott-state of one 
atom per site, the sampling procedure also represents all 
the many-particle correlations of the system. 

The joint probability distribution P(fi, . . . , r^v) for the 
positions of the iV atoms is given by the absolute square 
of the normalized many-particle wavefunction 

P(fi,...,r^) = |vl/(fi,...,rw)|2 (BIO) 

F(fi,...,f^)c.^ ^ |(/.,,(ri)p...|^,„(r^)P 

(Bll) 

Here the coordinates of the N atoms are denoted by 
(f 1, . . . , r^v). The summation runs over all possible A^- 
tuples (ji , . . . , jjv) of the state labels j = (1, . . . , refer- 
ring to the lattice sites j with precisely one atom per site. 
The Wannier wavefunctions for each lattice site (l>j in 
Eq. (|B11[) are assumed to have a negligible overlap with 
the neighboring sites, so that 0j(r)0fe(r) ~ 0, whenever 
j ^ k. Since we have exactly one atom per site, we may 
rewrite the joint probability distribution of Eq. (jBlip as 

P(fi,...,fAr)c.^ ^ |0l(f,jr..|0w(r,;„)|2 
ii...ijv 

(B12) 

where the summation runs over all possible permutations 
of the atomic coordinates i = (1, . . . , iV). Since there are 
iV! such permutations, we can express the joint proba- 
bility distribution in terms of the positions of an atom 
in the jth lattice site r^, instead of the position coor- 
dinates of the jth atom fj. Each A^! term contributes 
equally in the sum and wc obtain for the joint proba- 
bility distribution of the position of an atom in the lat- 
tice sites j = (1, . . . , A^) in the coordinate representation 
(ri,...,rAr) 

P'(ri, . . . , tat) ~ |0i(ri)|2 . . . |0^(r^)|2 (B13) 



By exploiting interactions between the atoms, one can 
tailor the incident field so that it drives specific linear 
combinations of collective modes, providing a desired col- 
lective response. In the text, for example, we consider 
phase modulated driving that can excite an array of iso- 
lated atoms arranged in a checkerboard pattern on an 
optical lattice. A spatial light modulator is employed to 
produce an incident field approximately of the form 

E+(r,t) = ej,i;oe*(*=^-"*)e'^(^'«\ (CI) 

where 

ip(x, y) = (/?i„ax sin(Ka;) sin(K?/). (C2) 

A field profile of this form, however, contains evanescent 
plane wave components whose transverse wavevectors ex- 
cede the the carrier wave number k. We therefore approx- 
imate this phase modulated field through the truncated 
Fourier expansion 

where km,ri = rriHiex + nuBy -\- Qm,n^zt Qm,n = 
\J (ri/c)2 — K?-{rv? -|- n^), and era,n is the normalized pro- 
jection of the vector hy onto the plane perpendicular to 
km,n- We have truncated the expansion for values of m 
and n for which „ > 0. 

Because our goal was to produce a phase modulated 
driving of the meta-atoms, we choose Cm,n so as to re- 
produce a phase modulated driving, i.e. 

d • E+ = ^ a„,„(d • e™,„)e'('''"--"-"*) (C4) 
^d-ej^e^^^^'^^e'C^^-"*), (C5) 

The coefficients Cm.n are obtained from the discrete 
Fourier transform of the phase modulation exp(z(^(a;, y)) 
and dividing by e„i^„ • d. 
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